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The Casimir and van der Waals interaction between two dissimilar thick dielectric plates is 

(N 

reconsidered on the basis of thermal quantum field theory in Matsubara formulation. We briefly 
review two main derivations of the Lifshitz formula in the framework of thermal quantum field 
theory without use of the fluctuation-dissipation theorem. A set of special conditions is formulated 
' under which these derivations remain valid in the presence of dissipation. The low-temperature 

■ behavior of the Casimir and van der Waals interactions between dissimilar dielectrics is found 

o, 

analytically from the Lifshitz theory for both an idealized model of dilute dielectrics and for real 



in 
o 



3 



X 



dielectrics with finite static dielectric permittivities. The free energy, pressure and entropy of 



the Casimir and van der Waals interactions at low temperatures demonstrate the same universal 
dependence on the temperature as was previously discovered for ideal metals. The entropy vanishes 

a, 

when temperature goes to zero proving the validity of the Nernst heat theorem. This solves the 



long-standing problem on the consistency of the Lifshitz theory with thermodynamics in the case of 
dielectric plates. The obtained asymptotic expressions are compared with numerical computations 
for both dissimilar and similar real dielectrics and found to be in excellent agreement. The role 
of the zero-frequency term in Matsubara sum is investigated in the case of dielectric plates. It is 
shown that the inclusion of conductivity in the model of dielectric response leads to the violation 
of the Nernst heat theorem. The applications of this result to the topical problems of noncontact 
atomic friction and the Casimir interaction between real metals are discussed. 
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I. INTRODUCTION 



Both the Casimir and van der Waals interactions are the quantum phenomena caused by 
fluctuating electromagnetic fields. In the framework of quantum field theory, these interac- 
tions can be described through the alteration by the material boundaries of the zero-point 
electromagnetic energy that pervades all of space. Thus, the Casimir and van der Waals 
forces are found to be closely connected with the fundamental properties of quantum vacuum 
(see the original Ref. [ij and monographs 2IO, U, 

In the last few years the Casimir interaction has been actively investigated in connec- 
tion with topical applications in extra-dimensional physics (where it provides an effective 
mechanism for spontaneous compactification of extra spatial dimensions) and in the bag 



model of hadrons 



a a. 



Several measurements of the Casimir and van der Waals forces 



between metal macrobodies were performed with increased precision using modern labo- 



ratory techniques 



s y, y, y y y 



13, 



14. 



15 



y y 



lo|. The experimental results 



were used to obtain stronger constraints on the Yukawa-type long-range interaction, pre- 



dicted bv 
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many extensions to 
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;he Standard Model, in the micrometer interaction range 
3Q. Consent with the dental application, in 



elementary particle physics, the use of the Casimir and van der Waals forces in nanotech 



nology 



26 



271 l28j ] . in quantum reflection and Bose-Einstein condensation 



noncontact atomic friction 



31 



32 



29 



30] 



and in 



was begun. 

The numerous applications of the Casimir and van der Waals forces and the extensive 
experimental work impose severe demands on the accuracy of theoretical predictions. The 



basic theory of these 
and his collaborators 



orces acting between real materials was developed by E. M. Lifshitz 
It expresses the free energy and force acting between two 



34, 



35 



macrobodies in terms of their frequency-dependent dielectric permittivities. The original 
derivation of the Lifshitz formulas for the free energy and pressure between two thick di- 
electric plates is based on the concept of a fluctuating electromagnetic field and uses the 
fluctuation-dissipation theorem. If the temperature of plates is not equal to zero, this con- 
cept includes both zero-point oscillations and thermal fluctuations of the field. In general, 
the Lifshitz derivation is applicable to both transparent and absorbing media. At present 
time, there are derivations of the Lifshitz formula based on the thermal quantum field the- 
ory in the Matsubara formulation (see Ref. p] for a review). They do not employ the 
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fluctuation-dissipation theorem but lead to results identical to those obtained with the use 
of this theorem. Thermal quantum field theory is immediately applicable to only transparent 
bodies, but, under certain conditions, dissipation can also be included. 

Quite independently from the Lifshitz formula, the Matsubara thermal quantum field 



theory was applied 



37 



38j for the investigation of the Casimir force acting between two 



plane parallel plates made of ideal metal (at zero temperature this case was considered 
already in Ref. Q). This is an exactly solvable problem of field quantization with the 
Dirichlet boundary conditions on the surface of the plates for the tangential component 
of electric field and with the identification condition in the Euclidean time variable. The 
obtained analytic expressions for the Casimir pressure, free energy and entropy were found 
to be in perfect agreement with thermodynamics. In particular, the magnitudes of the free 
energy and pressure were found to be monotonously increasing functions with the increase of 
the temperature. Entropy is nonnegative and goes to zero when temperature vanishes; this 
is in accordance with the third law of thermodynamics (the Nernst heat theorem) 2s|]. For 
the first time, the Lifshitz formula was considered to be in contradiction with the quantum 
field theoretical approach to the case of ideal metals since it leads to a different result in 
the limit of infinitely large dielectric permittivity. The situation was clarified by the so 
called Schwinger's prescription P, 40] whereby the limit of infinite dielectric permittivity 
in the Lifshitz formula must be taken prior to putting the Matsubara frequency equal to 
zero. Once that prescription is followed, the results obtained for ideal metals from the 
Lifshitz formula match those obtained from thermal quantum field theory with the Dirichlet 
boundary conditions. 

Currently the Lifshitz formula leads to problems when applied to the case of two plates 



41 



42, 



43, 



441 ] metallic plates were 



made of real metals at nonzero temperature. In Refs. 
characterized with the dielectric permittivity of the Drude model. In the limit of ideal metals 
this approach comes into conflict with the values of the Casimir free energy and pressure 
obtained using thermal quantum field theory with the Dirichlet boundary conditions. As 
was also proved in Refs. (2, 13 , the substitution of the Drude dielectric function in the 
Lifshitz formula leads to a violation of the third law of thermodynamics in the case of perfect 
crystal lattice with no defects and impurities (as discussed above, the quantum field theory 
approach obeys this law in the case of ideal metals). 

Another way to describe the realistic properties of a metal is to use the dielectric per- 
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mittivity of the free electron plasma model 



40 



47 



48( or to impose the surface-impedance 



boundary conditions |46j, |49j . In the case of ideal metals, both these approaches are in agree- 
ment with the results obtained from using thermal quantum field theory. They also obey 
the third law of thermodynamics. Recent experiments permit to test different theoretical 
approaches to the thermal Casimir force. It is significant that the first modern measuremen ; 
of the Casimir force Q| was found 5(3,13 ^° be in disagreement with Refs. 41. [l! 43. 44 1 



and consistent with Ref . 0] . The most precise and accurate experiments of Refs. |l 



exclude the theoretical approach of Refs. |4J 



tent with Refs. 



• 42,0Q 



42, 



43, 



16, 



441 ] at 99% confidence and are consis- 
. (Note that there is some controversy in the literature on the 



agreement of different approaches with thermodynamics and experiment which is reflected 
in recent Refs. 0,0].) 

As is evident from the foregoing, there are many questions in the Lifshitz theory of the 
Casimir and van der Waals interactions that remain to be answered. On the one hand, 
the consistency of the Lifshitz theory with thermal quantum field theory is not completely 
understood, and, on the other hand, the agreement of the Lifshitz theory with thermody- 
namics is called into question. An important point is that some thermodynamic aspects 
of the Lifshitz theory are still unknown even in the case of two dielectric plates for which 
consensus in the literature is achieved. One of the major problems realized 50 years back, 
but not resolved up to the present, is the elucidation of the low-temperature behavior of 
the Casimir (van der Waals) free energy, pressure and entropy between two dielectric plates. 
Without a resolution of this fundamental issue it would be impossible to settle the more 
complicated problems related to real metals. Even the correct way of practical computations 
for comparison of experimental results with theory would be uncertain. 

In the present paper we reconsider the thermal van der Waals and Casimir interac- 
tions between two thick dissimilar dielectric plates. We start with a brief analysis of two 
derivations of the Lifshitz formula in the framework of thermal quantum field theory in the 
Matsubara formulation. Special attention is paid to the conditions under which the quan- 
tum field theoretical derivations are applicable not only to transparent media but also to 
media with dissipation. We next consider the Casimir effect between two plates made of 
dissimilar dilute dielectrics with constant dielectric permittivities e^ k ' (k = 1, 2). In this 
case it is possible to develop the perturbation theory in the small parameters e^*' — 1 and to 
obtain the explicit analytic expressions for the Casimir free energy and pressure which are 
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exact in the separation between the plates and their temperature. The simple asymptotic 
behavior of the free energy, pressure and entropy at both low and high temperatures (short 
and large separations, respectively) is also found. The leading terms of the Casimir entropy 
at low temperature are shown to be of the second and third power in temperature, thus, 
demonstrating the agreement of the Lifshitz formula with thermodynamics in the case of 
dilute dielectrics. 

As a next step, we derive analytically the low-temperature behavior of the Casimir and 
van der Waals forces between dissimilar dielectrics using the frequency-dependent dielec- 
tric permittivities. For this purpose, perturbation theory in a small parameter which is 
proportional to the product of the separation distance between the plates and their tem- 
perature is developed. It is proved that the thermal corrections to the van der Waals and 
Casimir energy and pressure have an universal form and their asymptotic behavior can be 
calculated in terms of the static dielectric permittivities of the plate materials. Once again, 
the leading contributions to the entropy are shown to be of the second and third power in 
temprature which proves the agreement of the Lifshitz formula with thermodynamics in the 
case of dielectric plates with finite static dielectric permittivities. In the limit - 1 < 1 
the above results for dilute dielectrics are again obtained. Results for the special case of 
similar dielectrics are also provided. This means that the long-standing problem on the 
low-temperature properties of the Lifshitz theory in the case of dielectric plates is resolved. 

In what follows we perform the comparison between the obtained analytic results and 
numerical computations for dielectric plates made of silicon and vitreous silica. The cases of 
both dissimilar and similar plates are considered. The dielectric permittivities of silicon and 
vitreous silica along the imaginary frequency axis are found by means of the dispersion rela- 
tion using the tabulated optical data for both materials. In all cases the excellent agreement 
between the analytic and numerical results is observed below some definite temperature 
(separation) values. 

Finally we discuss the role of the zero-frequency term in the Matsubara sum for dielectrics. 
To calculate the van der Waals friction, Refs. [3^ . Q use the Lifshitz-type formula and 
describe the dielectric permittivity of the dielectric substrate by means of the Drude model 
with appropriately low conductivity. Formally, this leads to infinitely high permittivity of a 
dielectric at zero frequency and to the modification of the zero-frequency term of the Lifshitz 
formula. We demonstrate that such a modification results in the violation of the Nernst heat 
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theorem. Hence, it follows that the dc conductivity of dielectrics is irrelevant to the van der 
Waals and Casimir forces and must not be included in the model of dielectric response. This 
conclusion leads to important consequences for both the problem of atomic friction and for 
the Casimir and van der Waals interactions between real metals. 

The paper is organized as follows. In Sec. II, two main derivations of the Lifshitz formula 
for the free energy in the framework of thermal quantum field theory are briefly discussed. 
Special attention is paid to the restrictions imposed on the type of relaxation. Sec. Ill is 
devoted to the case of dissimilar dilute dielectrics. In Sec. IV, we derive the low-temperature 
behavior of the van der Waals and Casimir interactions between dissimilar dielectrics with 
frequency-dependent dielectric permittivities. Here, the consistency of the Lifshitz theory 
with thermodynamic requirements is proved. In Sec. V the obtained analytic results are 
compared with the results of numerical computations performed for some real dielectrics. 
Sec. VI demonstrates that the inclusion of the conductivity of dielectrics in the model of 
dielectric response leads to the violation of the Nernst heat theorem. In Sec. VII we present 
our conclusions and discussion touching on the topical problems of the noncontact atomic 
friction and the Casimir interaction between real metals. Appendices A, B and C contain 
the mathematical proofs of some statements used in Sees. IV and VI. 

II. DERIVATIONS OF THE LIFSHITZ FORMULA IN MATSUBARA QUAN- 
TUM FIELD THEORY 

We consider two thick dielectric plates (semispaces) with the frequency- dependent dielec- 
tric permittivities £^(w) and e^ 2 '(u), restricted by the two parallel planes z = ±a/2 with 
separation a between them, in thermal equilibrium at temperature T. 



The original Lifshitz derivation of his formula 34] was based on the assumption that the 
dielectric materials are characterized by randomly fluctuating sources of long wavelength 
electromagnetic fields. This concept includes not only the thermal fluctuations but also 
the zero-point oscillations of the field ■ Lifshitz based his derivation on the fluctuation- 
dissipation theorem using the properties of electromagnetic fluctuations. 

There are two derivations of the Lifshitz formula at nonzero temperature in the framework 
of thermal quantum field theory without use of the fluctuation-dissipation theorem. The 
first of them is due to M. Bordag (f| and is based on the scattering approach 5^. An 
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electromagnetic wave e lkz coming from the negative part of z-axis in the dielectric with the 
permittivity e^(u) will be scattered on the empty gap between the two semispaces and 
there occurs a transmitted and a reflected wave, 



(pi(z) ~ e ikz + si 2 e ikz , 



z— >— oo 



(pi(z) ~ sue ikz . (1) 

z— >oo 

A similar situation holds for the linear independent wave e~ tkz coming from the positive 
part of z-axis 

z^ — oo 

y> 2 (z) ~ e~ ikz + s 2 ie ikz . (2) 

z^oo 

The matrix {%} composed of the coefficients in Eqs. (JTJ) and (J2J) is unitary. 

For the evaluation of the free energy one has to use the Euclidean version of the field 
theory obtained by Wick rotation xq — > ix^ with the electromagnetic field periodic in the 
Euclidean time variable with a period (3 = h/(kBT), where ks is the Boltzmann constant. 
Starting from the usual representation of the free energy T in thermal quantum field theory 
at the one loop level, we arrive at the result (see Ref. p for details) 

1 „ \ r°° k ± dk ± 



2tt 



x 



lnsji(z6, k±) +lns^(i^, k±) . (3) 



Here sj\ are the elements of the scattering matrix for the two independent polarizations of 
electromagnetic field, k± = |kj_| is the magnitude of the wave vector in the plane of plates, 
and £t = 2nkBTl/h are the Matsubara frequencies. 

To obtain the explicit expressions for the elements of the scattering matrix, we assume 
that the electric field satisfy the Maxwell equations and the corresponding boundary condi- 
tions on the boundary surfaces z = ±a/2 (i.e., that the normal components of = e^'E^ 
and tangential components of are continuous). By solving the scattering problem with 
the demand that for infinitely remote plates both the free energy and pressure are equal to 
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zero, one obtains 



where 



s!i(«&>fci 



ft 



J. 



^ (1) ft + k?> e?> qi + k 



(!) J 2 l 



I 2 ) 



-2a(j; 



- ft k\*> -ft _ 



(2) 



2aq l 



(2) 



f 2 



(4) 



(5) 



cf V er 

An important point is that the Maxwell equations supplemented by the corresponding 

boundary conditions lead to a definite solution of the scattering problem, given by Eqs. (@J 
and only for I > 1. For I = 0, however, g = &o = an d, as a result, the respective 
system of linear algebraic equations has infinitely many solutions, i.e., it is satisfied for any 
s ii(0, k±) 5(|. The definite value of sj\(0, k±) is obtained by the use of the unitarity con- 
dition resulting in (s^O, k±)\ = 1 and dispersion relation leading to s^O, k±) = 1. We will 
return to this point below when discussing the role of dissipation. 

Substituting Eq. (jlj in Eq. (jSJ), we arrive at the Lifshitz formula for the free energy of 
the van der Waals and Casimir interaction 



k B T 
2vr 



1=0 



i - - 2 s u 



k i dk i 



x 



In 
+ ln 



'l-r^(^k x )rf\^k ± )e- 2 ^ 
l-ri 1} te,fc ± )rf ^ h k ± )e- 2 ^' 



(6) 



where the reflection coefficients are defined by 



(fe) 



(fe) , 7 W ' 

e\ >qi + k\ 



(6, h 



(fe) 



ft 



*f } + ft 



(7) 



The second field theoretical approach to the derivation of the Lifshitz formula at nonzero 
temperature without recourse to the fluctuation-dissipation theorem is based on the direct 



summation of the free energies of all photon oscillator modes |49| (at zero temperature this 
m e th „d was Pr0 p„ 8 e d ta Ref, HQ; see als0 the ions in Ref, QflH). 



Equations for the determination of the frequencies of oscillator modes between two dis- 
similar plates are obtained from the Maxwell equations with the corresponding boundary 
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conditions on the boundary planes z = ±a/2. They are given by 6Ct 161] 



An (u>, k 



o)q\ ( "" = 0, 



A ± (u,k ± ) = (kW+q) (k ( V+q) e a(1 



-aq 



where 

The solutions of Eq. 
mode is given by 



to 



ki-sW(u;)-, q 



k, - -7T. 



(8) 



(9) 



cr V cr 

can be denoted as u\ n and n . The free energy of one oscillator 



■7*l,n( a . T ) 



+ k B Tln \l-e~ 



fcsTln 2 sinh 



k± ,n 



2k B T 

After summation over all quantum numbers, we obtain 

k B T 



(10) 



J 7 (a, T) 



2tt 



E 



k i dk i 



In 2 sinh 



k i ,n 



2k B T 



+ In 2 sinh 



k i ,n 



2k B T 



(11) 



Eq. pi)) can be identically represented by the use of the argument theorem in terms of 
the quantities A\\(u,k±) and A_i_(u,k±) from Eq. (jHJ) having their roots at the frequencies 
of the oscillator modes. After transformation of the branch points into poles by means of 
integration by parts and calculation of the residues at all imaginary Matsubara frequencies 
we rewrite Eq. (jllj) in the form (see Ref. H for details) 

k B T " 



2tt 



1=0 



E ( i - 5*- 



k i dk i 



[lnA||(^,fc±) + lnA ± (z^,fcj 



(12) 



Expression (fT2*J) is infinite. To remove the divergences, we subtract from the right-hand side 
of Eq. (fT2*j) the free energy in the case of infinitely remote plates (a oo). Then Eq. © 
with the notation (J2J) is reobtained. 



Both quantum field theoretical derivations of the Lifshitz formula discussed above are 
immediately applicable for dielectrics described by real dielectric permittivities eW(w). In 
this case the unitarity condition in the first derivation is valid and the photon oscillator 
modes in the second derivation are real quantities. Both derivations, however, under certain 
conditions, can be generalized to the case of media with dissipation Q, 62]. In so doing a 
medium under the influence of electromagnetic oscillations can be represented as a set of 
oscillators 

d ^ + ^ + i °° x{t) = m (13) 
where 7 is a damping parameter and f(t) = exp(-iu)t) is the harmonically changing 
external force. Eq. (JTHJ) has the solution x(t) = x^exp^—iut) where i w = x{ u )fui an d 

xH = -5-^ 2 ( 14 ) 

is the susceptibility of the system. The latter is connected with the dielectric permittivity 
of a medium 

e(u) = i + \(u), (15) 

71 

where g is the oscillator strength. 

The dielectric permittivity (j!5|) is complex and describes an absorption band (or bands if 
one considers several oscillators with different parameters) of finite width and amplitude at 
nonzero frequency. Note that in thermal equilibrium net heat losses are absent on average 



62 



63j |. One can conclude that for such media all processes of absorption are balanced by 



the respective processes of emission, and the global unitarity condition remains valid. 

Regarding the second field theoretical approach described above, the oscillator frequency, 
determined from the equation l/x(w) = 0, becomes complex when 7 7^ 0. In this case the 
free energy is not given by Eq. (jlljl (which is already clear from the complexity of the right- 
hand side of this equation). For the complex dielectric permittivity given by Eq. (fTKj) the 
correct expression for the free energy is obtained from the auxiliary electrodynamic problem 



and leads to Eqs. (fEZJl and © (see Ref. 62] for details). The qualitative reason for the 



validity of this statement is that the free energy depends only on the values of the dielectric 
permittivity along the imaginary frequency axis [i.e., only on e(i£)] which are always real. 
Note that these considerations are not applicable to real metals described by the Drude 
dielectric function because in this case the proper frequency uj of the respective oscillator 
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in Eq. (fT3|) turns into zero and the absorption band is shifted to zero frequency and achieves 
an infinite amplitude. As a result, dissipation leads to heating of a metal and the unitarity 
condition is violated (see Refs. |49l l56l|). 



III. CASIMIR EFFECT FOR TWO PLATES MADE OF DISSIMILAR DILUTE 
DIELECTRICS 



In this section we anticipate that the dielectric permittivities of both plates (k = 1, 2) are 
constant and equal to = 1 + ijk, where rjk <C 1. The assumption that the dielectric per- 
mittivity does not depend on frequency implies that the separation distance a is rather large 
(in fact larger than the characteristic absorption wavelength) and, consequently, relativistic 
effects are essential aa |. Thus, in this case we are dealing with the Casimir effect. 

For convenience in analytic and numerical calculations, we introduce the dimensionless 
variables ( and y given by 

6 2o£, 







rl, y = 2q l a, 



(16) 



where £ c = c/(2a) is the so called characteristic frequency, r = Aixk,BaT/ (he), and qi was 
defined in Eq. (jSJ). Then the Lifshitz formula © takes the form 



F{a,T) 



her 

327r 2 a 3 



E 



ydy 



x jln 
+ In 



i=o x / " £ 

r^\Ci,y)rf\( l ,y)e-y 



(17) 



r { t\tuyVl\tuy) 



(2), 



Here the reflection coefficients (|IJ) are expressed in terms of new variables ()16|) as follows: 



4 k) y-\/y 2 + (?(z 



ri k \d,y) 



2 (J® 



rf(0,y) 



y 2 + C!(e\ k) -i)-y 



y' + Q{ef ] -l)+y 



(18) 



Note that ef^ = e^{i$j) = e^ k \i(iC, c ) . In the case of constant dielectric permittivities, which 
are under consideration in this section, e, = e^ k \ i.e., being independent on I. 
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Substitution of eW = 1 + rjk in Eq. (J 17)) and expansion of the logarithms in powers of the 
small parameters rji, r]2 (preserving all powers up to order 3 inclusive) leads to 



where 



2567r 2 a 3 



1=0 



X 



(19) 



" Ci (2 + 20 + I - f ) + C 2 (2 - f ) Ei(-O), 



^ (2y 6 - Cfy 4 - Cfy 2 + Cf) 



e [ + 4 + 12 + 24 24 



+ < i+ !-S) ei (- 



-co, 



and Ei(;?) is the exponential integral function. 

Summation in Eq. (fTT?|) can be performed in the following way: 



fi(r) 



E (1 - bio) uci) 



(1 + r)e T - 1 rV [(1 - r)e 2r - 4re r - r - 1] 



(e- - l) 2 

00 

+I>i(co. 



2(e r - l) 4 



/ 2 (r) =£(1- ^0) h 



(0) 



(20) 



(21) 



(22) 



(1 + r)e T - 1 r 2 e T [(5r - 3)e 2r + 20re r + 5r + 3] 
+ (e- - l) 2 + Y2{? - l) 4 

tV [(r - l)e 4r + 2(13r - 5)e 3r + 66re 2r + 2(13r + 5)e T + r + 1] 



24(e r - l) e 



+ E F2 ^), 



1=1 



12 



where 



./•i.(\/) = sf(2-f )Kii-,„ 



(23) 



2 24, 

Thus, the exact (in r) expression for the Casimir free energy for the configuration of two 



^(G) = C? ( 1 + ? -!r) Ei (-0)- 



dissimilar dilute dielectric plates is given by 

T{a,T) - 



(24) 



256?r 2 a 3 

where /i )2 (t) are defined in Eqs. lj2Tl Hj23|) . 

It is not difficult to find the asymptotic behavior of Eq. ()24j) at r -C 1. For this purpose 
the remaining sums on the right-hand side of Eqs. (|21|l . (|22} are calculated with the help of 
the Abel- Plana formula ^, [" 



°° / I \ POO 

Y.y 1 - 2 5i °) = I F{t)dt 

F(it) - F(-it) 



poo 

+ i dt~ 
Jo 



(25) 



where F(z) is an analytic function in the right half-plane. Taking into account that ^1,2(0) 
0, these sums are equal to 



00 poo 

yX 2 (C0 = / F la {rt)dt 

1=1 J ° 

F 1>2 {irt) - F li2 {-irt) 



POO 

+ i dt- 
Jo 



2irt 



(26) 



Using Eq. (j23j) one obtains 



Fi(rt)dt 



16 
15? 



F 2 (rt)dt 



128 



(27) 



'0 ±UI Jo 

The second integral on the right-hand side of Eq. 
From Eq. fl23) it follows: 



105r 

can be calculated perturbatively. 



F x (iTt) - F±(-iTt) = -2inrH 2 + AirH 3 + 0(r 4 ) 
F 2 (irt) - F 2 (-irt) = -itr¥ + 2irH 3 + 0(r 4 ). 



(28) 
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Substituting Eqs. (j2*Tj) and (|2H|) into Eq. (|2T)j). we arrive at 

16 C(3) 



2=1 

oo 



15r 
128 



(29) 



i=i 



105r 4vr 2 



T 

120 



0(r 4 



where £(2) is the Riemann zeta-f unction. After expansion of the remaining terms on the 
right-hand sides of Eqs. (j2*Tj) . (J2~2~j) in powers of r, one obtains 

46 C(3) o 7r 3 



/ 2 (r) 



15r 
338 



360 



+ 0(r 4 



(30) 



r 



105r ' 4vr 2 ' ' 360 
The substitution of Eq. (|3*U|) in Eq. (|2*ljl leads to the final result 



F{a,T) 



hcr\\T]2 46 
~256vr 2 a 3 [l5 
tji + t/a /338 , C(3) 



C(3) r 3 
2tt 2 



+ 



III 

360 
f 0(r 5 



(31) 



105 ' 4tt 2 ' ' 360 / 

The asymptotic behavior of the Casimir free energy at r ^> 1 (high temperatures or large 
separations) also can be obtained from Eq. (|24|). With the proviso that t>1, Eqs. (|2*Tj) 
and (|2*2*Jl lead to /i(r) = f%{r) = 1 up to exponentially small corrections. As a result 



Vl +V2 



(32) 



The same expression follows from the zero-frequency contribution to the free energy in 
Eq. (HZ!) 
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V1V2 



(2 + Vl )(2 + m ) 



(33) 



after integration with respect to y and expansion in powers of rjx and r] 2 . Note that at zero 
frequency only contributes to Eq. (jHBJ), because, according to Eq. (fTHJ) . r±\o,y) = 0. 

Now let us consider the Casimir pressure and entropy for the configuration of two dis- 
similar dielectric plates. The Lifshitz formula for the Casimir pressure presented in terms of 
dimensionless variables (JTfij) has the form 



P(a,T) 



her 



X 



32vr 2 a 4 

1=0 

rP(( h y)rf\( h y) 



ey - yK'(Cu y) e» - r£%, y)rHG, v) 



,( 2 ) 



— Sin 

2 



+ 



y 2 dy 
ri 1) (G,l/)r? ) (G,») 



(34) 



14 



Expanding in powers of small parameters 771, r\2 by the same procedure as in the case of 
the free energy, we obtain the expression exact in r 



P(a,T) 



ftCT VlV2 
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f \ V1+V2 , s 
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where Pi,2( T ) are given by 
Pi (t) = 2 + 4 
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In the limiting case r< 1, Eqs. (|33jl . (|3^|) lead to 

^071772 
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This expression can also be obtained as 

P(a,T) 



dTja, T) 
da 



(38) 



where !F(a,T) is given by Eq. (J3TJ). 

At high temperatures (large separations) it holds r ^> 1 and from Eqs. (J33j) . (pjj) it 
follows: 



P(a,T) 
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32yra 3 

Using the above procedure, one can obtain the Casimir entropy for the configuration of 
two dissimilar dilute dielectric plates. It can be found also as 

&F(a, T) 



S(a,T) 



dT 



where !F(a,T) is given by Eq. (|24|1. In the limit r < 1 the result is 



5(a,T) = 



3fc B C(3)r7ir/ 2 r 2 
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As is seen from Eq. (J41J) . S(a, 0) = 0, i.e., the Nernst heat theorem is satisfied. 

The obtained results (p?T]) . (jHZj), and (JIT]) are used as the tests in the next section where 
the low-temperature behavior the Casimir and van der Waals interactions between plates 
with arbitrary (not dilute) dielectric permittivities is studied. 

IV. LOW-TEMPERATURE BEHAVIOR OF THE VAN DER WAALS AND 
CASIMIR INTERACTIONS BETWEEN DISSIMILAR DIELECTRICS: ANALYT- 
ICAL RESULTS 

We will now look at the configuration of two parallel plates made of real dissimilar di- 
electrics described by frequency-dependent dielectric permittivities. A distinguishing feature 
of dielectrics is the finite value of their static dielectric permittivity Eq = £^(0). It is com- 
mon for dielectrics that the dielectric permittivity is practically equal to its static value 
within some frequency region [0,0^]. At higher frequencies e^ k \u) is smaller than . 

We start with the Lifshitz formula (fTTj) . (fTHj) which is applicable to any dielectric plates 
at any separation (note that separations a should be large enough so one can ignore the 
atomic structure of the plates, i.e., a should be larger than 2 or 3nm). 

Applying the Abel-Plana formula (|25)l . we can identically rearrange Eq. (fTTj) to the form 

F(a,T)=E(a)+AJ r (a,T). (42) 

Here E(a) is the energy of the van der Waals or Casimir interaction at zero temperature 
given by 

he 

~ 3l 

where we introduce the notation 



E{a) = 32^J d( J c dyf &> y )> (43) 



f((,y) = y\\n I - r\' {U)r\\U)e- y 



+ ln 



l-r { l\C,y)r { l\c,y)e- y ]}. (44) 



The second contribution on the right-hand side in Eq. ()42j) is the thermal correction to the 
energy 

A _ m , iher f 00 F(irt) - F(-irt) 
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where the function F(x) is defined by 

F(x) = j dyf(x,y). (46) 



The behavior of the energy ()43|) at both short separations (the limit of the van der 
Waals forces) and large separations (the limit of the Casimir forces) is well-studied [36]. 
At short separations (a <C A ) the dependence of the dielectric permittivities on frequency 
is important and the asymptotic behavior of the energy is given by E(a) = —H/(12ira 2 ), 
where H is the Hamaker constant JfJ [jf3|. At large separations (a ^> Ao) £^(u;) can be 
approximated by their static values Eq . As a result, the energy takes the form E(a) = 
— $?(eq~\ e^)/a 3 , where \& is some tabulated function of two variables ((I l^. 

Importantly, the low-temperature behavior of the thermal correction ()45|) under the con- 
dition r < 1 has a universal form valid at both short and large separations. The qualitative 
explanation for this fact, proven below in detail, is that for dielectrics at sufficiently low 
temperatures the Matsubara frequencies contributing to the thermal correction belong to 
the region where practically do not depend on the frequency and are equal to their static 
values Eq . 

To prove the universal low-temperature character of the thermal correction, we use the 
Ninham-Parsegian representation for the dielectric permittivities along the imaginary fre- 
quency axis 120, E3, 

Here, the parameters Cj are the absorption strengths for different dielectrics (k = 1, 2) 
satisfying the condition, 

£cf = 4 fc) -l, (48) 

j 

and uij k ^ are the characteristic absorption frequencies. Although Eq. (}4"Tj) is the approximate 
one, it gives a very accurate description for many dielectrics [65]. 

To obtain the asymptotic behavior of the thermal correction to the energy given by 
Eq. (|4"H|) at r < 1, we substitute Eq. (jUJ) in Eq. (jUJ) and expand the function f(x,y) in 
powers of x = rt. The subsequent integration of this expansion with respect to y from x to 
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infinity (see Appendix A) leads to 

(1) , (2) , 2 (1) (2) , (1) _ l)( (2) _ ^ 

(4 1) + 1)(4 2) + 1) 2(e«+4 2) ) 
-iax 3 + 0(x 4 ). (49) 

On the right-hand side of Eq. (J49)) we separated the third-order contribution with a real 
coefficient a. It cannot be determined at this stage of our calculations because all terms 
in the expansion of f(x,y) in powers of x contribute to this coefficient [the value of a is 
found below in Eq. ([6ip]. As is seen from Eq. (}4*9*|) . only the static dielectric permittivities 
contribute to F(ix) — F(—ix) in leading order and, consequently, to the thermal correction. 
Substituting Eq. (|4l?j) in Eqs. (jj3J) and and performing the integration, we arrive at 



the result 



fxc 

F{a,T) = E{a) - -— (50) 



32vr 2 a 3 

c(3)(4 1) + 4 2) + 24 1) 4 2) ) ( £ «-i)(sf-i; 



(4 1} + i)(4 2) + i) 8^(4^ +4 2 ); 



r 3 - C 4 r 4 + 0(r 5 ) 



where C 4 = a/240. 

So far we have considered the free energy. The pressure can be obtained, using Eqs. (J55J) 
and (JoT)j). as 

P(o, T) = P (a) - — - [C 4 r 4 + 0(r 5 )] , (51) 

where Po(a) = —dE(a)/da is the pressure at zero temperature. Our aim is to determine the 
value of C4. To attain this, the pressure is expressed directly by the Lifshitz formula (|34jl . 
Applying the Abel-Plana formula ()25j) in Eq. (|34|). we represent the pressure in the form 

P(a,T)=P (a) + AP( ffl ,T), (52) 

where the thermal correction to the pressure is 

and the function $(x) = $i|(x) + <f>±(x) is given by 

$ , ^ _ r™ y 2d y r {^( x >y) r \i±( x >y) ^ 

' Jx e y — AP ± (x, y)r^p ± (x, y) 



18 



First, let us determine the leading term of the expansion of in powers of x. By 

introducing the new variable v = y/x one arrives at 



$,(*)=* 3 / t^l^±A^ , (55 ) 
J i e vx — r ± (x, v )r ± (x, v ) 

Note that the reflection coefficients r±\x,v) depend on x only through the frequency de- 
pendence of according to Eq. (JJ7J). Expanding in powers of x in Eq. (f55|). we obtain 



3 /• oo i; 2 dT;ri 1) (0,T;)ri 2) (0,T;) 4 , 
$±(x)=x i / ^ - - ^ - - + 0(x 4 ), 



1 1 -r ( l\0,vyi>(0,v) 



(2), 



(56) 



where, according to Eq. (|T8*j) . 



V 2 + £q^ — 1 — V 

v 2 + 4 fc) - 1 + U 



(57) 
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£ o £ o ~ £ o £ o \ x° 



Integration of Eq. with account of Eq. (J57jl results in 

$ ± (x) = 
From Eq. (|55|) it follows 



fc o + V fc o 



+ 0(x 4 ). 



(58) 



$ i (zz) - $ i f-zx) = — i 1 



(1) (2) ./_(1)_(2) _(1)_(2) , 3 

e o + £ o +y £ o £ o ~ £ o £ o \ x 6 
£ o ^ V fc o 



+ 0(s 5 



(59) 



The expansion of 3>ii(x) from Eq. ()54|) in powers of x is somewhat more cumbersome. It 
is presented in detail in Appendix B leading to the result 



$ii (ix) - $\\(—ix 



-z< 1 + 



£ o + V £ o "* ) ( £ ^ + £ ^ 
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Adding Eqs. and (pjUj) and integrating the obtained result according to Eq. ([53]). we 
find the coefficient C4 in the thermal correction to the energy ()50|) and pressure ()51|) 
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This expression gets much more simplified for the case of two plates made of one and the 



same dielectric (e 



(i) 



-(2) 



£ 0)- 



C, 



720 



Thus, for two similar dielectrics the free energy takes the form 



F{a,T) = E{a) 



he 
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C(3)( g p-1) 2 
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(63) 



720 



(y/EE - 1) { £ l + ^V^o - 2) r 4 + 0(r 5 



Another limiting case is given by dilute dielectrics. Putting Eq = l+r]k in Eqs. (I5U|) . (foT|) 
and (|61|) and expanding in powers of the small parameters rjx, r] 2 , one arrives at precisely 
the same thermal corrections (terms proportional to r 3 and r 4 ) as were obtained for dilute 
dielectrics in Eqs. (f3~Tj) and (J3ZJ). Note that Eqs. (f3~Tj) and (J37)) also contain terms independent 
of r having a physical meaning of the energy and pressure at zero temperature. These terms 
coincide with E(a) and Po(a) from Eqs. ()50j) and (|51|1 only in the relativistic limit of large 
separations. This is because Eqs. (}3*Tj) and (|37j) were obtained under the condition that the 
dielectric permittivities are constant. However, as indicated above, the thermal corrections 
obtained in such a way have a universal character and are valid at any separation with the 
constraint r< 1. 

Eqs. flpTIj) . (JBTj) and ()63|) solve the vital issue about the thermodynamic consistency of the 
Lifshitz formula for the case of two dielectric plates. Using Eqs. (|5U|l and (jjDj) , the entropy 
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of the interaction between plates takes the form 



S(a,T) 
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(64) 



3C(3)(e ~~ ± A c o 

where the coefficient C4 is determined in Eq. (jfiljl . As is seen from Eq. ()64j) . the entropy 
vanishes when the temperature goes to zero as it must be in accordance with the third law 
of thermodynamics (the Nernst heat theorem). 

In the limiting case of two similar dielectrics Eq. can be rearranged as follows: 



S(a,T) 
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(65) 
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In the limiting case of two dilute dielectrics the expansion of Eq. ()64)1 in powers of rji and 
1] 2 coincides with Eq. (I41j) . 

As is seen from Eqs. fl64|) and ()65j) . in the limit r — > (T — ► 0) the lower order contri- 
butions to the entropy are of the second and third powers in the small parameter r. In this 
manner for dielectrics at low temperatures the entropy obeys the same universal law which 

66 1 and real Q| metals. Recall that ideal 



was previously found for ideal 



38, 



39 



45 



metals (i.e., two plates with the Dirichlet boundar y co nditions E 4 = 0) lead to an ex actly 

38]. As was shown in Refs. 39l |45| . 
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solvable model in Matsubara quantum field theory 
in the case of plates made of ideal metal the entropy at low temperatures is given by 

2tt 2 



3k B ((3) 2 



1 - 



-r + 0(r 2 ) 



(66) 



135C(3) 

This demonstrates that the thermal quantum field theory approach is in perfect agreement 
with thermodynamics for ideal metals. It should be pointed out that the expansion coef- 
ficients of the free energy, pressure and entropy in powers of r in the case of ideal metals 
cannot be obtained as a straightforward limit — > 00 in Eqs. ()5(J|) . ([HT]) and ()64j). The 
mathematical reason is that it is impermissible to interchange the limits r — >• and e — > 00 
in the power expansion of functions depending on e as a parameter. 

In the case of large separations (high temperatures) r 3> 1 and the approximation of 
static dielectric permittivities is applicable. With this condition the major contribution is 
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given by the zero-frequency term of the Lifshitz formula (|17|) 
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(the other terms being exponentially small). Integration in Eq. (JBTj) leads to 
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where Li n (z) is the polylogarithm function. In a similar manner at r ^> 1 the pressure is 
given by 

k B T. 



P(a,T) 



rLi a 
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(69) 



In the limiting case of dilute dielectrics Eqs. 
respectively. 



and (jHHj) coincide with Eqs. (J32j) and 



V. COMPARISON BETWEEN ANALYTIC AND NUMERICAL RESULTS 

In this section we compare our analytic results in Eqs. (J5U|) and (J51|) for the thermal 
corrections to the free energy and pressure with computations performed for both dissimilar 
and similar dielectrics using the complete tabulated optical data and the Lifshitz formulas 
f!17|) . (jHlJ). This enables us to illustrate the applicability regions of the above asymptotic 
expressions for different plate materials and to gain an impression on the temperature de- 
pendence of the free energy and pressure within a wide range of parameters. As an example 
we consider two dielectrics (Si and vitreous SiC^) which differ in their dielectric properties. 

The dielectric permittivity of real dielectrics along the imaginary frequency axis can be 
obtained through the dispersion relation 

= 1 + - / duj y > (70) 
^ Jo u + f 

Here Ime(w) = 2ra 1 (u;)n 2 (co'), where rii(u) = Ren(w), U2(uj) = Imn(u), and n(ou) is the 
complex refractive index tabulated in Ref. js2|. For the dielectric Si (of high resistivity 
p = 1000 Q cm) and vitreous Si0 2 the dielectric permittivities were computed by the use of 
Eq. (|7nj) in Ref. [68| . The obtained results are presented in Fig. 1 for Si (line 1) and for SiC>2 
(line 2). The flat steps in both lines should be extended for all frequencies < £ < 10 10 rad/s 
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(the lowest frequency indicated in Fig. 1). As a result, the static dielectric permittivities 
are equal to £q = 11.66 (for Si) and = 3.84 (for SiC^). 

In Fig. 2, we compare the computational results for the thermal corrections to the energy 
AjF(a,T) = jF(a,T) — E(a) (a) and pressure AP(a,T) = P(a,T) — Po(a) (b) in the con- 
figuration of two plates one made of Si and another one of SiC>2 at separation a = 400 nm 
as a function of the temperature. Short-dashed lines are obtained by the use of the Lifshitz 
formulas (JT7J), (J4*3~)) with the static dielectric permittivities £q . Solid lines are com- 

puted using the same Lifshitz formulas with the complete frequency-dependent dielectric 
permittivities of Si and SiC>2 presented in Fig. 1. Long-dashed lines show our asymptotic 
expressions for A J 7 (a, T) and AP(a,T) on the right-hand sides of Eqs. (|qT7|) and (joTj) . re- 
spectively, with a coefficient C 4 defined in Eq. (JBTj) . 

As is seen from Fig. 2a, the thermal correction to the energy can be calculated with 
the static dielectric permittivities of Si and Si0 2 at temperatures below 100 K. At higher 
temperatures the solid line computed by the use of the tabulated optical data departs from 
the short-dashed line computed with the static dielectric permittivities. At T < 60 K the 
asymptotic expression for AT{a 1 T) in Eq. (|5()j) with the coefficient (|6*Tj) leads to the same 
results as the original Lifshitz formula. Quite similar situation takes place for the thermal 
correction to the pressure (see Fig. 2b). The numerical results obtained by the use of the 
Lifshitz formula with static coincide with those obtained using the tabulated optical 
data at T < 100 K (solid and short-dashed lines). As to the asymptotic expression for the 
thermal correction to the pressure [Eqs. (|51|) and (jfilj) ]. it becomes applicable at T < 50 K. 

To give a more comprehensive idea on the application regions of the obtained asymptotic 
expressions, in Fig. 3 we present the same information, as in Fig. 2, for the case of two 
similar plates made of vitreous SiC>2 at separation a = 450 nm. As before, in Fig. 3a the 
thermal correction to the energy is plotted and in Fig. 3b to the pressure versus temperature. 
As is seen from Fig. 3a, all three theoretical descriptions of the free energy (in terms of the 
tabulated data, static dielectric permittivity and the asymptotic one) become applicable at 
T < 85 K. If we consider the thermal correction to the pressure (Fig. 3b), different models 
of £ are applicable at T < 100 K whereas the asymptotic expression ()51|) with the coefficient 
(162}) coincides with them at T < 65 K. The distinguishing feature of Fig. 3b is the intersection 
of the solid and long-dashed lines. It is explained by a greater deviation of the solid line, 
computed by the use of tabulated data, from the short-dashed line, computed using static 
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because of the presence of a second flat step in the frequency dependence of e^ 2 \i^) 
(see line 2 in Fig. 1). 

In Fig. 4 we present the dependences of the thermal correction to the energy (a) and 
pressure (b) for two Si plates at separation a = 300 nm computed in the framework of the 
same approaches. For two Si plates within the separation and temperature region under 
consideration the lines computed with frequency- dependent and static almost coincide 
(see Figs, 4a,b). The asymptotic expressions in Eqs. (|U3*|) and ()51|) with the coefficient (|fi2"|) 
are applicable at T < 45 K and T < 35 K, respectively. In Fig. 4a the asymptotic expression 
for the thermal correction to the free energy achieves its maximal value at about T w 105 K. 
This is explained by the relatively large static dielectric permittivity of Si. 

Comparing all the above figures, one can conclude that the obtained asymptotic expres- 
sions for the free energy (containing two terms of order r 3 and r 4 ) have a wider application 
range than the asymptotic expressions for the pressure containing only one term of order 
t 4 . At the same time, by decreasing the separation distance one can widen the range of 
temperatures where our asymptotic expressions are applicable. The remarkable feature of 
Figs. 2-4 is the monotonous increase of the magnitude of all thermal corrections with the 
increase of temperature for real dielectrics (see solid lines overlapping with our asymptotic 
expressions in the applicability region of the latter). This confirms the same conclusion 
made in Refs. |69| based on qualitative thermodynamical considerations. 

VI. ROLE OF THE ZERO-FREQUENCY TERM IN THE MATSUBARA SUM 
FOR DIELECTRICS 

According to Eq. (|17jl . the free energy of the van der Waals and Casimir interaction 
is represented by the Matsubara sum from zero to infinity. The zero-frequency term in 
Eq. (H7J, 

j rp poo 

F (a,T) = -^ ydy In l-r^O, y)rj 2) (0,y)e-" , (71) 

(k) 

is of special interest. For dielectrics with finite static dielectric permittivities e from 
Eq. (JIEJ) it follows: 



,(*) 

+ I 



r^) = ^_-^ rl k \0,y)=0. (72) 
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Eqs. (fTTj) and (|72jl were already used in Eqs. (f3l?j) and (fBTj) to obtain the asymptotic expres- 
sions for the free energy at high temperatures where the contributions from all terms in the 
Matsubara sum with I > 1 are exponentially small. 

Note that the second equality in Eqs. (|72j) is somewhat analogous with the same equality 



in the case of real metals described by the Drude dielectric function [4l( . In reality, however, 
for metals the second reflection coefficient at zero frequency is r\\(0,y) = 1, i.e., equal to 
its physical value. Together with r±(0,y) = 0, this leads to the violation of the Nernst 



heat theorem for perfect crystal lattices with no impurities [46] and to contradictions with 
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experiment |l5l.ll6l.ll7l 3,Q|. On the contrary, for dielectrics ri k \o,y) in Eq. (|72j) is larger 



than its physical value at normal incidence [the latther is equal to (y — l)/(y £ o^ + -Q] 
and, as was demonstrated in Sec. IV, the Lifshitz formula incorporating Eqs. (fT2*j) is in 
perfect agreement with the Nernst heat theorem. 

We now turn to a problem of outstanding importance which arises when one includes the 
dc conductivity of dielectric materials into the model of their dielectric res pon se. What this 
means is that, instead of dielectric permittivities = e^(i^i), one uses p^L lo^ 

e f » . = + ^f- = ef + SUM. (73) 

Here Oq is the dc conductivity of the plate materials and (3^ = 2ha^ j '{ksT). It is 
common knowledge [70| that the conductivity of dielectrics depends on the temperature as 
cr^ ~ exp(— /T) where is determined by the energy gap A^) which is different for 
different materials. It cannot be too highly stressed that for dielectrics the conductivity 
at constant current is very low. To take an example 7l| , for SiC>2 at T = 300 K it holds 
p( 2 ) 1CT 12 and, thus, negligible for all I > 1. From physical considerations the inclusion 
of the term (5^ k \T)/l into the model of the dielectric response (f7H|) seems of dubious value 
since (5^ k \T) quickly decreases with decrease of T and, thus, remains negligible at any T. 

In spite of this, the substitution of Eq. (|75j) into the Lifshitz formula (fTTj) leads to Eq. (fTTj) 
with different value of one of the reflection coefficients at zero frequency than in Eq. (|72p: 

rf\0,y) = l, rf\0,y) = 0. (74) 

The change from Eq. (|72j) to Eq. (j74j) has far-reaching consequences. To analyze them, we 
calculate the free energy J- {a, T) with the dielectric permittivities e z identically rearranging 
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Eq. (JT7J) to the form 
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where the reflection coefficients f^\_ are obtained from Eq. (|18|) by replacing e\ K> with e 

Now we expand the last, third term on the right-hand side of Eq. (J75J) in powers of the 
small parameters (3^/1. Combining the zero-order contribution in this expansion with the 
second term on the right-hand side of Eq. (J7SJ), one obtains the free energy J 7 (a, T) calculated 
with the dielectric permittivities ej . Calculating the first integral on the right-hand side 
of Eq. ()75|). we arrive at 



T(a, T) = J 7 (a, T) 
k B T 



(76) 



167ra 2 



C(3) - Li 3 



(2) 



(4 1} + i)(4 z; + i) 

In this formula R(a,T) is of order 0(P^/l) and it stands for the first and higher-order 
contributions in the expansion of the third term on the right-hand side of Eq. (|75|) in powers 
of (3^/1. The explicit expression for R(a, T) is given in Eqs. (jClj) and ()C2)) of Appendix C. 
As is shown in Appendix C, R(a, T) exponentially goes to zero with the decrease of T. 

Eq. (J76|) leads to the important conclusion about the thermodynamic inconsistency of the 
Lifshitz formula for dielectrics if one includes the dc conductivity in the model of dielectric 
response. Substituting Eq. (|75jl into Eq. (jjO]), we obtain the entropy for the plates with the 
dielectric permittivities (|75jl in the form 



-(2) 



+ R(a,T). 



S(a, T) = S(a, T) 
k B 



(77) 



167ra 2 



C(3) - Li 3 



4 1} -i)(4 z; -i) 



(2) 



4 1} + i)(4 ZJ + i 



(2) 



dR(a,T) 
df 
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where S(a,T) is the entropy for the plates with the dielectric permittivities given by 
Eq. (JES). 



In the limit T 



S(a,Q) 



it follows: 

k B 



167ra 2 



C(3) - Li 3 



(4 i; -i)(4 2) -i) 



-(1) 
-(1) 



4 i; + i)(4 2) + i) 



> o. 



(78) 



Eq. (|78|) depends on the parameter of the system under consideration (the separation dis- 
tance a) and implies the violation of the Nernst heat theorem. Thus, the dc conductivity 
of a dielectric is irrelevant to the origin of the van der Waals and Casimir forces and must 
not be included in the models of the dielectric response. The neglect of this rule results 
in the violation of thermodynamics. Physically it is amply clear that for high-frequency 
phenomena like the van der Waals and Casimir forces the behavior of dielectric materials 
at low frequencies is described by the static dielectric permittivities (the approach used by 
E. M. Lifshitz and his collaborators [35 1). The results of this section provide the necessary 
theoretical background for this conclusion. 



VII. CONCLUSIONS AND DISCUSSION 



In the foregoing, we have reconsidered the Lifshitz theory of the Casimir and van der 
Waals interaction between dielectrics on the basis of thermal quantum field theory in Mat- 
subara formulation. As was shown above, there are field theoretical derivations of the Lifshitz 
formula without use of the fluctuation-dissipation theorem which permit generalization for 
the presence of dissipation. The special analysis demonstrates, however, that only a partic- 
ular dissipation is compatible with the Lifshitz formula, those without net heat losses and 
with balanced processes of absorption and emission. 

In this paper, we have analytically solved the long-standing problem of the low- 
temperature (short separation) behavior of the thermal corrections to the Casimir energy 
and pressure between dissimilar dielectric plates and demonstrated that the Casimir (van 
der Waals) entropy vanishes when the temperature goes to zero. This was done, first, using 
the idealized model of dilute dielectrics, and, then, for real dielectrics with finite static di- 
electric permittivities. The free energy, pressure and entropy of both the van der Waals and 
Casimir interactions between dielectrics demonstrate at low temperatures the same universal 
temperature dependence which was previously discovered for ideal metals. This proves the 
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thermodynamic consistency of the Lifshitz theory in application to dielectric plates. The 
obtained asymptotic expressions were compared with the results of numerical computations 
for both similar and dissimilar real dielectrics (silicon and vitreous silica), and were found 
to be in excellent agreement. 

Special attention was paid to the role of the zero-frequency term in Matsubara sum in 
the case of dielectric plates and to the extrapolation of the dielectric permittivity along 
the imaginary frequency axis to zero frequency. As was proved above, the inclusion of the 
conductivity at constant current in the model of dielectric response leads to a modification 
of the zero-frequency term of the Lifshitz formula and to the violation of the Nernst heat 
theorem. The conclusion was made that the conductivity of dielectrics at constant current 
(low but formally nonzero at nonzero temperature) is irrelevant to physical phenomena 
described by the Lifshitz theory. This conclusion leads to far-reaching consequences for 
both the problem of noncontact atomic friction and for Casimir interaction between real 
metals at nonzero temperature. 

As to the problem of atomic friction, the observed effect is many orders of magnitude 
larger [3l| than the van der Waals friction between metals caused by the vacuum fluctuations 

n 

and thermal photons (in the experiment of Ref. |31( the gold-coated tip and substrate were 
used). In Refs. 32|, |54j it has been proposed that measurements performed on metal films 
are strongly affected by the underlying dielectric substrate. To calculate the friction due to 
the substrate, Refs. 0, S] use the Lifshitz-type formula including the low dc conductivity 



in the model of dielectric response. Then, the large value of friction coefficient is obtained 

m «B 

proposition of Refs. |32|, |54| would not be in agreement with the Nernst heat theorem. It is 
also qualitatively clear that the low-frequency behavior of the dielectric permittivity in the 
region below ~ 600rad/s cannot cause the large friction effect at characteristic frequencies 
from 3.75 x 10 14 rad/s to 1.5 x 10 17 rad/s (in the experiment of Ref.jsj] separations vary 
from 1 to 400 nm). 

Regarding the application of the obtained results to real metals, the case of two dielectric 
plates with including the dc conductivity turns out to be analogous to two metal plates 
described by the Drude dielectric function. In both cases the region of low frequencies 
contributes significantly to the description of high-frequency phenomena resulting in con- 
tradictions with thermodynamics (see Introduction). It is well known that in the region of 
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the normal skin effect, described by the Drude dielectric function, dissipation leads to the 
heating of a metal (see Ref. for details). For dissipation of this kind, as discussed above, 
the field theoretical derivation of the Lifshitz formula is inapplicable. In this connection 
it is not surprising also that the results obtained from the Lifshitz formula combined with 
the Drude model are in conflict with the conclusions obtained for ideal metals using the 
thermal quantum field theory approach with the Dirichlet boundary conditions. The pro- 
posed measurement of the Casimir force between both metals and semiconductors at large 
I — n 

separations [72] is aimed to resolve this contradiction experimentally. At large separations 
(high temperatures) the classical limit of the Casimir effect is achieved where forces acting 



between real and ideal metals practically coincide 
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74j . In our opinion, this important 



experiment could bring the final verification of the results obtained from the thermal quan- 



tum field theory approach (see also Refs. [75. [76J] where two short-separation experiments 
are proposed with the same aim). 

It is worthwhile to note that in real experiments the most frequent configuration is not 
the two parallel plates but a large sphere above a plate. There is a prediction in recent 
literature [z^f that for ideal metals in the case of open geometries the thermal correction 
to the Casimir force is not as suppressed as in the parallel plate geometry in the limit of 
zero temperature. By this reason it is of much interest to derive on fundamental grounds 
the low-temperature behavior of the Casimir free energy between real bodies (metallic or 
dielectric) for any open geometry. 

To conclude, the above reconsideration of the Casimir and van der Waals interactions 
between dielectric semispaces not only confirms the mutual agreement between the thermal 
quantum field theory approach, Lifshitz formula and thermodynamics, but also suggests 
ways for the resolution of other complicated problems in modern applications of quantum 
electrodynamics. 
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APPENDIX A 



In this Appendix we prove Eq. ()49|) in Sec. IV. Let us present Eq. (|44jl in the form 



f((,y) = f\\((,y) + f±((,y), 



(Al) 



where 



(A2) 



f l± (C,y) = 2/ ^ [l - r«(C,y)rj 2 l(C,y) e ^ . 

It is easy to check that f±(C,y) does not contribute to the first term on the right-hand side 
of Eq. (|49p. Substituting Eq. (|4*7jl in /||(C, 2/) and expanding in powers of x = rt one obtains 



J 1 ) i J 2 ) , 



1 — Tq Vq e 



+ 



fc ' ^ fc o fc 



r r e 



+ l) (e^ + l)y(l-r^e-y 



x 



(A3) 



+ 2 



1 



Sc 



+ 



1 



v 4 2) + ij ^ 

where the static values of the reflection coefficients are found from Eq. (fTSj) 



.(2)' 



r r ^ x 2 + 0(x3) 



1 - rJ 1) rS 2) e-» 



4 fe) + 1 



< i. 



(A4) 



Note that for simplicity we consider only one oscillator in Eq. (|47|) and omit index j. The 
obtained results, however, are valid for any number of oscillators. 

As a next step, we should integrate Eq. (jA3|) according to Eq. (|46|). Integration of the 
first term on the right-hand side of Eq. (jA3|) results in 



roc 

Zi(x) = / ydyln ( 1 — ^Vq e~ 9 



n=l 



'1 + nx)e 



n" 



(A5) 



Expanding (|A5|) in powers of x and summing the obtained series, we arrive at 

Z 1 (x) = -Lis (rfVf ) 

- y In (l - r^rf) +0(x 3 ). (A6) 
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Integration of the second term on the right-hand side of Eq. (jA3|) contains the integral 

OO (1) (2) -y 

2 / j„ '0 '0 e 



Z 2 (x) = x dy 



y yl — r^Vo^e ^ 

oo 

2 (rfV^E^-nx). (A7) 



-x 

n=l 



Integration of the third term on the right-hand side of Eq. (|A3j) leads to the integral 

2 / J„.„. '0 e 



poo 

Z 3 (x) = x 2 dyy 

J X 



1 _ T ,( 1 ) r ( 2 ) y 

i r r e 



E(w) n (i+ :f - . (as) 



n=l 

Expanding in powers of x in Eq. (jA8j) and performing the summation, we obtain 

4 (1) (2) 

Z 3 (x) = x 2 U 2 (rfVf > ) - X - -^ g (2) + 0(z 5 ). (A9) 

1 r r 

Let us now calculate the different contributions to the quantity (|4"9*|) . According to 
Eq. (JMJ) 

^i(ia?) - Zi(-ix) = 0(x 3 ). (A10) 



From Eq. ()A7|) it follows: 

(1) (2) 

Z 2 (zx) - Z 2 (-ix) = tirx 2 T ° (2) + 0(x 3 ). (All) 
Finally, Eq. flM]) leads to 

Z 3 (za;) - Z 3 (-ix) = 0(x 5 ). (A12) 

As is seen from Eqs. (IA10|) - (jA12|) . only Eq. ()A11|) contributes to the leading term on the 
right-hand side of Eq. (|49p. As to the frequency dependence of the dielectric permittivity 
in accordance to Eq. (|47jh it contributes only to the fifth order term [see the last term on 
the right-hand side of Eq. (jSHJ) and Eq. (fA"l2| ]. 

Assigning the numerical coefficient to Z 2 as in the second term on the right-hand side of 
Eq. (HHJ) and using Eq. IjAllJ) . we arrive at the result 

(i) , (2) , 2 £ {1) e {2) mr {1) r {2) 

f ^ - -h*> - if:^ + <v>- (a 13) 

In terms of a notation (jA4|) . this coincides with Eq. ()49|l. 
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APPENDIX B 



This Appendix is devoted to the derivation of Eq. (|6U|) in Sec. IV where the function $n (ix) 
was used in the calculation of pressure. As was shown in Appendix A, the dependence of 
the dielectric permittivity on frequency contributes to the expansion of F(ix) in powers of 
x, used in the calculation of the free energy, starting from only the 5th power [i.e., to <&\\(ix) 
starting from the 4th power in x}. Here we are looking for the lowest (third) order expansion 
term of $>\\(ix). Because of this, it is possible to disregard the frequency dependence of e 
and describe the plate materials by their static dielectric permittivities. 

Let us identically rearrange Eq. (J54j) by subtracting and adding the two first expansion 
terms of the function under the integral in powers of x 

r,| (x, y)rf (x, y)e~ y 



$11 fx) 



dy 



r^\x 1 y)ru Z> (x,y)e y 
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dy 



1 _ r ( 1 ) r ( 2 ) p -'i / 

i r r e » 



2 ' 



We consider the first integral on the right-hand side of Eq. (jBlj) written in terms of a 
new variable v = y/x 

r\{\0,v)rf\0,v)e-™ o r^r^e'"" 
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where, in accordance with Eq. (|18jh 



fB2l 
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(fc) 

4 ^ 



^ 2 + 4 fe) -i 
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(B3) 
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The leading expansion order of Qi(x) from Eq. (|B2|) in powers of x is 

rS 1) (0,'u)rf | 2) (0,'u) 



l-rP(0,v)r\? } (0,v) 
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2 ' ' 
1 — 7"g 7g 
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(4 iJ + l)(4' ; + l) (l-rfrf)' 
The explicit calculation of the integral in Eq. (jB4|) leads to 
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Now we calculate the second integral on the right-hand side of Eq. (IBljl . 



Q 2 (x) 



Tq Tq e v 
1 - rj 1} rj 2) e-^ 



EK 



.(1)„(2) 




r, 



n (2 + 2m; + n 2 x 2 )e~'" : 



n=l 



Expanding Eq. ()B6|) in powers of x, one obtains 

Q 2 (x) = 2Li 3 (r«rf ) - X - - 



n- 



^ _(1)_(2) 
'0 '0 



_ r (l) r (2) 1 ( X )• 

The integral contained in the third term on the right-hand side of Eq. (jBl|) is 
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calculated 



Q 3 (x)=x 2 / dy - ^f" n2 (B8) 
1 - rj VJ J e-» 



2 ' ' 



1 — r^V^e 



It can be expanded in powers of follows: 

J1U2) J1W2) 

Q 3 (*) = x 2 -^L - s» r ° + 0(**). (B9) 
1 - r o V o M - rfVfj 

Finally, according to Eq. (|B1|) . the function under consideration is given by 

(4 1) + i)(4 2) + i: 

where Qi(x), Q2(x), and Qz(x) are found in Eqs. ()B5|) . ()B7|) . and ()B9|) . respectively. It is 
notable that the contributions of the third power from Qi(x) and Qz(x) in Eq. ()B10|) cancel 
the second and third contributions from Q\{x). 

Using Eq. (jBlOJ) . one arrives at Eq. (|6U|) after some identical rearrangements. 



$H (x) = Qx(x) + Q 2 (x) - ° n ° : J° Qs(*), (BIO) 



APPENDIX C 



The quantity R(a, T) was introduced in Eq. (|76|) of Sec. VI and has the following explicit 
form: 

R(a,T) = R w (a,T) + i? (2) (a,T) + \((3 {k) /l)' 



where 



i? (1) (a,T) 



dyy e 



k B T 

1 1 >/r + vc; j, -i; 



(2-e«)Cf-2y 2 
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.(2) 



r\\Cuy) 



(CI) 



(C2) 



v / 2/2 + C 2 (£ (i)_ 1) + y l-rf ) (C J ,y)ri a) (Ci,y)e-»_ 

and R( 2 \a,T) is obtained from R^(a,T) by interchanging of the upper indices (1) and (2). 
In this Appendix we demonstrate that R(a,T) vanishes exponentially when T — > 0. 
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Let us consider the integral with respect to y from Eq. (|C2|) . expand the integrated 
function in powers of r (we recall that (i = rl) and restrict ourselves by the main contribution 
to the result at r = 0: 



dy- 



ye 



,(2) 



(C3) 



where was defined in Eq. (JMJ). With account of Eq. (jSJ), Eq. J0H1) can be rearranged 
to the form 
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Substituting Eq. (d in Eq. (|H2j) . we find 
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It is readily seen that 

- In (1 - e" nT ) + 
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lnr + 1 - Inn + 0(r 2 



(C6) 



The substitution of the leading term on the right-hand side of Eq. (IC6J) in Eq. ()C5|) results 
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Taking into account that (3^' ~ {\/T) exp(—b^/T) (see Sec. VI), we get the conclusion that 
the temperature dependence of (a, T) is determined by the term 

i?«(a,T)~e- fe(1) / T lnT, (C8) 

i.e., both R^\a, T) and its derivative go to zero as exp(— U- 1 ' /T) when the temperature goes 
to zero. 

In perfect analogy to R^(a, T), the same conclusion is obtained for R^ 2 \a, T). The terms 
of the second and higher powers in 13^ [see Eq. (|UTJ)] go to zero even faster than R {k \a,T) 
when T — > 0. In this way, with account of Eq. (|C1J) . we have proven the exponentially fast 
vanishing of R(a, T) with decrease of the temperature. 
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FIG. 1: Dielectric permittivity of Si (line 1) and vitreous Si02 (line 2) along the imaginary 
frequency axis as a function of the logarithm of frequency. 
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FIG. 2: Magnitudes of the thermal corrections to the energy (a) and pressure (b) in configuration 
of two plates one made of Si and another one of SiO"2 at a separation a = 400 nm as a function of 
temperature calculated by the use of different approaches: by the Lifshitz formula and tabulated 
optical data (solid lines), by the Lifshitz formula and static dielectric permittivities (short-dashed 
lines), by the asymptotic expressions in Eqs. Q5U|) . (|51|) and (|61|) (long-dashed lines). 
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FIG. 3: Magnitudes of the thermal corrections to the energy (a) and pressure (b) in configuration 
of two plates made of vitreous Si02 at a separation a = 450 nm as a function of temperature 
calculated by the use of different approaches: by the Lifshitz formula and tabulated optical data 
(solid lines), by the Lifshitz formula and static dielectric permittivities (short-dashed lines), by the 
asymptotic expressions in Eqs. (|STj) . and (|p^)l (long-dashed lines). 
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FIG. 4: Magnitudes of the thermal corrections to the energy (a) and pressure (b) in configuration 
of two plates made of Si at a separation a = 300 nm as a function of temperature calculated by 
the use of different approaches: by the Lifshitz formula and tabulated optical data (solid lines), 
by the Lifshitz formula and static dielectric permittivities (short-dashed lines), by the asymptotic 
expressions in Eqs. (|51j). ()62j) and (|^3*)) (long-dashed lines). 
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